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Abstract 

In this paper, we study the dynamics of shearfree dissipative gravi- 
tational collapse in modified Gauss-Bonnet theory of gravity - so called 
f{G) gravity. The field equations of /(G) gravity are applied to shear- 
free spherical interior geometry of the dissipative star. We formulate 
the dynamical as well as transport equations and then couple them to 
investigate the process of collapse. We conclude that the gravitational 
force in this theory is much stronger as compared to general relativity 
which indicates increase in the rate of collapse. Further, the relation 
between the Weyl tensor and matter components is established. This 
shows that for constant f{G) the vanishing of the Weyl scalar leads 
to the homogeneity in energy density and vice versa. 

Keywords: Modified Gauss-Bonnet gravity; Dissipative fluid; Gravitational 
collapse. 



1 Introduction 

Recently, a lot of work has been done on alternative theories of gravity for 
the identification of dark energy (DE) which is said to be responsible for 
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accelerated expansion of the universe. These theories are attractive due to 
their consistency with the astrophysical observations and local gravitational 
experiments.^-* The simplest generalization to general relativity (GR) is f{R) 
gravity in which / is an arbitrary function of the Ricci scalar R. Although, it 
is the simplest modification of GR but it is not possible generally to formulate 
such f{R) model which is consistent with the solar system tests. Many 
conditions for the validity of f{R) models have been proposed. 

Another modified theory of gravity "the modified Gauss-Bonnet theory" 
was proposed by several authors. "^^^^^ In this theory, GR is modified by in- 
troducing some arbitrary function f{G) in Einstein-Hilbert action, where 
G = R^-4:R^"'R^^ + R^"''^^R^^^s is the Gauss-Bonnet invariant. This modifi- 
cation (consistent with observational constraints) is related to string-inspired 
dilaton theory. Cognola et al.^-* found that such model can produce transi- 
tion form matter dominated era to accelerated phase. 

The dynamics of the dissipative shearfree gravitational collapse is an im- 
portant issue. In GR, this problem was initially formulated many years 
ago by Misner and Sharp^^ for adiabatic collapse and by Misner^^^ for non- 
adiabatic collapse. The dissipative process during the collapse of massive 
star is important due to the energy loss during formation of neutron star 
or black hole.^^^ It is shown^^^'^^^ that gravitational collapse is a dissipative 
process, so the effects of the dissipation must be included in the study of col- 
lapse. Herrera and Santos^^-* discussed the dynamics of shearfree anisotropic 
radiating fluid. Chan^^^ studied the realistic model of radiating star with 
shear viscosity. 

Herrera^^) explored the inertia of heat and its role in the dynamics of dis- 
sipative collapse. Herrera et al.^^^ also derived the dynamical equations by in- 
cluding dissipation in the form of heat flow, radiation, shear and bulk viscos- 
ity and then coupled with causal transport equation. They investigated that 
shearfree condition for quasi-static (slowly evolving) non-dissipative systems 
in Newtonian limit leads to homologous collapse. Di Frisco et al.^^^ investi- 
gated the dynamics of spherically symmetric charged viscous non-adiabatic 
gravitational collapse. Sharif and his collaborators ^9)~25) j^avc discussed the 
dynamics of charged and uncharged viscous dissipative gravitational collapse 
in GR as well as in f{R) gravity. 

In this paper, we discuss the dynamics of the dissipative shearfree gravita- 
tional collapse in /(G) gravity. We use the field equations of /(G) gravity 
developed by using covariant gauge invariant (CGI) perturbation approach 
with 3-1-1 formalism. The collapsing matter is taken as dissipative isotropic 
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fluid bounded by the shearfree interior geometry. This spacetime is matched 
with the Vaidya exterior solution by using Darmois junction conditions 
The plan of the paper is as follows: In the next section, we describe the 
field equations in f{G) gravity and matching conditions. We formulate the 
dynamical as well as transport equations and then couple them in section 
3. In section 4, we discuss the relation between the Weyl tensor and matter 
variables. The last section provides discussion of the results. 



2 Field Equations in Modified Gauss-Bonnet 
Gravity 

The Einstein-Hilbert action for the modified Gauss-Bonnet gravity is 

5 = /Ay=5(^±P + £„), (1) 

where g is the determinant of the metric tensor g^j^^, R is the Ricci scalar, 
f{G) is an arbitrary function of Gauss Bonnet invariant G, Cm is the matter 
Lagrangian and k is coupling constant. The variation of this action with 
respect to the metric tensor gives the field equations 

G^u — nTfj,^ + -g^vf — "^FRR^^ + AFR^^R^^ — 2FR^^sujRl^'^ 

- AFR^^s.B?' + 2i?V^V,F - 2^^,V'F - AR^V^V^F 

- ARZV^V^F + m^.V'^F + Ag^^R^^V^VsF - 4i?^^,5 V^V^F, (2) 

where F — df{G)/dG. The energy-momentum tensor for dissipative fluid is 
T,j,u = {p + p)V^Vy - pg^u + q^Mu + QuV^,, (3) 

where p, p, and are density, pressure, four velocity and radial heat 
flux, respectively. 

In general, it is very difficult to handle the field equations. We use a 
simphfied form of these equations derived by Li et al.^^^ through CGI pertur- 
bation using 3-1-1 formalism and the effective energy-momentum approach. 
Since there are nonlinear terms in R^^ and R/^uXa which are expressed in terms 
of dynamical quantities, so the field equations and Gauss-Bonnet invariant 
are simplified as follows 

G,,^k{T,, + T;^,), (4) 
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where is the Gauss-Bonnet correction term. Using CGI approach, the 
components of are evaluated in terms of the dynamical quantities as^^^ 



= iQ(/-FG) + ^(fi-3M/)(Fe + V2F)), (5) 

-P^ = iQ(/-FG) + ^(fi-3vI/)F-^^](0F + V^F)^, (6) 

= ^-^(fi-3v&)(V,F-ieV,F) + ^F0C,). (7) 

Here, V is spatial covariant derivative, 6 = V^^ is expansion scalar and the 
quantities fl, are expressed in terms of dynamical quantities as follows 

n = -(e + ^e^-vM^), (8) 



^ = -^{(d + e' + R-WA^), (9) 

C = -^^ + VV^. + V'^c.M-, (10) 

where R is the Ricci scalar of 3D spatial spherical surface, cr^j, = V(^:j,) — 
- l^h^^ (where h^^ = g^^ - V^V^) is the shear tensor, A^, = V^.^V 
is four acceleration and uj^^ = V[^.^] + V[^V^] is the vorticity tensor. The 
Gauss-Bonnet invariant G in CGI approach is 

G' = 2(^i?2_^M-^^^)^ (11) 

where 

R = -20 - + 2V''Af, - R, 
R^^R,, =1{Q' + 00' + lo') + ^(0 + 0')^ - ^(0 + 1^')^'^,. 

Now we assume that the fluid ([3]) is bounded by the spherical boundary 
S. The metric interior to S is assumed to be co moving and shearfree in the 
following form 

ds^ = XHe - Y\dr^ + rMe^ + sin^ edcf"), (12) 
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where X and Y are functions of t and r. The corresponding four velocity 
and heat flux take the form 



V>' = X-H\, q^ = q6^i, (13) 

where V^qf^ = 0. The expansion scalar for the fluid sphere is = Using 
Eqs.dSD, ©, (113) and ([H]), the field equations yield 



8.iq + q-)XY^ = 2 | y - ^ - || ) , (15) 



7r(p + p^)r2 = (_) +ri_ + l^)+2 



'Y'V 2 /y X'\ X'Y' fY^ ^ 



Y ) r \ Y X } XY \X 

Y f y\^ XY 

n n , F" 1 /r X'\ X" 

8n{p + p^)rY' = r'{—-^yj + - + ^ j + X 



where p^, ]P and g*^ correspond to Gauss-Bonnet contribution to GR. Here 
dot and prime indicate derivative with respect to time and radial coordinate, 
respectively. Making use of Eqs.([5])-([7j), these turn out to be 

X (^-V^^h (18) 
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Kp^ = ^-{GF-f) 



1 1 /6F fY'^"^ 



8 / r X F 1 /j^// x''^ X' 





Y"\ •• 






Yr ^ 




fx" 




\'x' 


X + 



3 \ Xy Y XY X2r2 3^2 \^ X X rX 
^^^^ V^fV (19) 



\ xr 



4 I QY [Y'V AY' Y"\ lYF' 



3y2 ' \Y J Yr Y J \YX 

8FX f YY' YX' Y' \ 
^ YY\XY^^X^~^j' ^ ' 

The Gauss-Bonnet invariant is found from Eq. flTT]) as 
G 




It follows from Eq.(IT5ll that 



(21) 



8vr(g + g«)F^ = ^. (22) 

If q+q^ = 0, then 0' = and we have the condition for homogenous collapse. 
The Misner-Sharp mass becomes^-* 

-(M) = -^---— 1. (23) 
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We assume that the exterior region to boundary surface S is described by the 
Vaidya spacetime which represents the geometry of the radiating star given 
by 

rfs^ = (1 - ^^^)du^ + 2dudR - R\de^ + sin^ ed(j)'^). (24) 
R 

Applying the same procedure as discussed 

ini9)-25), we find that the conti- 
nuity of the first and second fundamental forms over the boundary surface 
E gives the continuity of the gravitational masses across E {i.e., m{r,t) = 
M(z/)) and p + p'^ =^ (g + q'^)Y. This implies that the effective pressure is 
balanced by the effective radial heat flux on the boundary. 



3 Dynamical and Transport Equations 

In order to study the dynamics of the field equations, we use the Misner- 
Sharp formalism^). Here velocity of the collapsing fluid is given by 

U = rDtV, (25) 

where U < for radially inward motion of the fluid and Dt = From 
the definition of the mass function fl25D. we can write 



(yy -fi^u^^'^y=E. (26) 



Y \ rY ^ 

This is known as the energy of the collapsing fluid. The proper time derivative 
of the mass function (1231) leads to 



..YYY l(Y\^ YY^X lYY'^ .Y'Y' Y' 
^*"^ = ^ '^ + 2lx) -^C^^2XY^-'yX-Vx^ 

' (27) 

Using Eqs.(Il5D, ([M]), ([2l]) and ([25]), this can be written as 

Dtm = -Att [{p + q^)U + (g + q^)YE\ {rYf (28) 

which give the variation of the fluid energy inside a sphere of radius rY . 
Since [/ < 0, so the first term on right side of this equation increases the 
energy of the system while the negative sign with second term implies the 
leaving of energy from the system. 
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We also define proper radial derivative Dr = where R = rY is the 
proper areal radius of the sphere. The proper radial derivative of Eg. (1231) 
leads to 

Y f .Y'Y" ,Y" fY'Y 3r^fY'^^ 
Y' ^X'Y^ r^Y'Y^ .Y'Y 3r^ f yV'. , , 

Using Eqs.(ITl]), (JH]), ([2l]) and ([25]) in ([29]), it follows that 

Dnm = An {{p + p^) + (g + q'')Y^ {rYf. (30) 

This equation describes the variation in energy between two adjacent layers 
of the fluid inside the spherical boundary. The Gauss-Bonnet term effects 
the density and heat flux. Since heat flux is directed outward due to U < 0, 
so the Gauss-Bonnet term causes the system to radiate away effectively. 

The acceleration of the infalling matter inside the spherical boundary is 
obtained by using Eqs.( IT6]) . ( 123]) and ( 125]) as 

W = -(^ + 4>r(p + ,°)(ry))+^^. (31) 

The r component of the energy-momentum conservation, (T^^ +T'-^^^);/3 = 0, 
implies that 

^{p + g^)' + ^{p + p"" + P + q"") + {q + f)^ + 5(g + q"")^ = 0- (32) 

Using the value of ^ from this equation in Eq. (l3T]) with ([23]) and ([25]) . we 
obtain the dynamical equation 

{p + p° +P + q^)DtU = -{p + p^ +P + q^) [m + 47r(p + q^)R^] 

- [Dr{p + g^)] - E[hY{q + q'')^ + YDt{q + g^)]. (33) 

This equation has the "Newtonian" form, i.e., force = mass density x 
acceleration, where mass density = p + p'~' + p + . This implies that 
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Gauss-Bonnet term effects tlie mass density due to fiigfier curvature. Tfie 
first square brackets on riglit side is tlie gravitational force, wliosc Newtonian 
part is m and relativistic part is p + q'-^ , lience tlie Gauss-Bonnet term effects 
gravitational force. The second term is the hydrodynamical force, it resists 
against collapse because -Di?(p -I- q*^) < 0. The last square brackets give the 
contribution of heat flux to the dynamics of the collapsing system. Here the 
first term is positive {U < 0, {q + q^) > 0), implying that outflow of heat 
flux reduces the rate of collapse by producing radiations in the exterior of 
the collapsing sphere. The effects of Dt{q + q'^) will be explained below by 
introducing the heat transport equation. 

Now we discuss the transportation of heat during the shearfree collapse of 
radiative fluid in f{G) gravity. For this purpose, we use the equation derived 
from Muller-Israel-Stewart phenomenological theory for non-adiabatic fluids 
28),29) jg known that Maxwell- Fourier law^°^ for the heat flux leads to 
diffusion equation which indicates perturbation at very high speed. For the 
relativistic non-adiabatic fluids, this is given by Eckart^^^ and Landau. To 
resolve this problem, many relativistic theories^^^"^''^ have been proposed. 
The common point of all these theories is their validity for non-vanishing 
relaxation time. These theories provide the heat transport equation which is 
hyperbolic equation. 

The transport equation for the heat flux in this case reads^^^ 

rh^-V%,x + r - ->^h^''{Tu + TA,) + ]^KT^ (^-^^ r, (34) 

where /i'*'^ is the projection tensor, q = q + q*^, t is relaxation time, T is 
temperature and k, is thermal conductivity. For the interior spacetime, this 
equation reduces to 

T^[{q + q'')Y] + iq + q'')XY' 



-KiTAy - 



KT\q + qG)Y^ r . 
2 ^kT'^' 



- 3 



TYY{q + q^ 



(35) 
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Using Eqs.( 12^ and (15U]) . this implies 

YDt{q + q'') = -KT 



G^_ ..rr^DtU kT Y{q + q^) ^ , rU 



tE tY t \ rY 

T 

dt \kT^ 
3UY{q + q^) 



(36) 



2R 

In order to see the effects of heat flux on the dynamics of collapsing 
sphere in the modified Gauss-Bonnet gravity, we couple the dynamical and 
heat transport equations given by Eqs. (l33|) and (!36|) . This coupling leads to 

EkT' 

(p + + p + - a)DtU = Fgravil - a) + Fhyd + 
E{q + q^)Y AE{q + q^)YU nT^q + q^)Y d ( t \ 



T 



R 2XE dt \kT^ 







Here 



. EI<l^. (37, 



Fgrav = -{p + +p + q^){m + A7i{p + q^)R^) ^ 



hyd 

a 



R^' 

E'[Dnip + q'')], 
kT 



t{p + p^ +P + q^Y 

Now we analyze the effects of the Gauss-Bonnet term on the dynamics of 
the collapsing radiative fluid. Since the presence of Gauss-Bonnet term in a 
(thermal coefficient) effects the value of a and the factor (1 — a) being the 
the multiple of Fgrav would effects the value of Fg^av In the definition of a, 
K, r and T (as mentioned after Eq. flSll) ) are positive for a real dissipative 
gravitating source ^"^^ and {p + p) > 0, for collapsing fluid to satisfy the weak 
energy condition. For p*^ and q'^, we have two cases 

• When p*^ > and g*^ > 0, then value of a in this case will be less as 
compared to GR and value of (1 — a) will be larger (as compared to 
GR) and consequently gravitational force will be stronger as compared 
to GR and rate of collapse will be increased. For example if f{G) = /o, 
a positive constant, Eqs. f|T8|) -f l20|) yield p'^ = —p^ = y and g*^ = 
and we have p*^ = y > and rate of collapse will be increased. 
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• When p'^ < and g*-^ < 0, then value of a in this case will be larger as 
compared to GR and value of (1 — a) will be less (as compared to GR) 
and consequently gravitational force will be weaker as compared to GR 
and rate of collapse will be decreased. For example if f{G) = /o < 
in above case, then we have p"-^ = y < and rate of collapse will be 
decreased. This can be verified for the more general f{G) models. 

If a approaches to a critical value of 1 then the inertial and gravitational 
forces of the system become zero and the system will be in hydrostatic equi- 
librium. Notice that a appears in the system due to its thermal conductivity 
and temperature. The value of temperature for which a — )■ 1 is equivalent 
to the expected amount of temperature that might be reached during the 
radiative collapse in supernova explosion Thus a increases during the 
supernova explosion if the relaxation time r 7^ 0. When a crosses the criti- 
cal value, the gravitational force plays the role of antigravity and expansion 
would occur in the system. 



4 Relation Between Weyl Tensor and Matter 
Variables 

In this section, we establish a relation between the Weyl tensor and density in- 
homogeneity which helps to extract some information about the gravitational 
arrow of time. The Weyl tensor leads to the Weyl scalar = C^"^^'^C^i.ao- 
which is given by 

= 7^ - 2R''I^R^p + ^R^ (38) 

where TZ = R^'^^" R^^Xa is the Kretchman scalar. Using the field equations 
in Eq.( H6|l in appendix, Eq.( l38|) yields 

e = m-^R\p + p''), (39) 

where 

^ -.R\ (40) 



V48 

From Eqs.dSO]) and it follows that 



Due = AtiR^ 



{q + f)y^-DR{p + p'')^ 



(41) 
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For /(G) = /o, a constant, Eqs.([IH])-([2n]) yield = -jP = f and = 
0, which imphes that for constant f{G), we have ACDM, model. Assuming 
that the fluid is non-dissipative, Eq.( BT]) yields 

Due + ^R'Dnip + ^) = 0. (42) 

If we take D^e = Dn (^-j^R^^ = 0, using regular axis condition, i.e., R ^ 0, 

we have C = 0. Thus Eq. fj42|) implies that D^ip = 0. This means that the 
conformal flatness condition produces homogeneity in the energy density and 
vice versa. Equation f H2|) gives the relationship between the energy density 
inhomogeneity and the Weyl tenor for perfect fluid. This is an important 
expression for the Penrose proposal about the gravitational arrow of time^^-* . 
The tidal forces associated with the Weyl tensor make the fluid more inho- 
mogeneous as the evolution occurs, indicating the sense of time. We would 
like to mention here that such a relationship is no longer valid if f{G) is not a 
constant, like local anisotropy of the pressure, dissipation and electric charge. 
In this case, we see from Eq. (l4ip how f{G) affects the link between the Weyl 
tensor and energy density inhomogeneity, suggesting that f{G) should enter 
into the definition of gravitational arrow of time. 



5 Summary 

In this paper, we have considered a new framework in which an arbitrary 
function f{G) of Gauss-Bonnet invariant is introduced in the Einstein-Hilbert 
action to account for the dynamics of the shearfree dissipative gravitational 
collapse. The corresponding field equations are used to study the shearfree 
dissipative interior geometry of the star. The matching of the interior to the 
exterior Vaidya geometry implies that the effective pressure is balanced by 
the effective radial heat flux on the boundary. Using the Misner deflnition, 
we have formulated the dynamical equation ( 133|) which shows that the Gauss- 
Bonnet correction increases (decreases) the inertial mass leading to increase 
(decrease) the rate of collapse. 

To discuss the transportation of heat during the shearfree collapse of 
radiative fluid, we have formulated heat transport equation. Also, we have 
coupled the dynamical equation with the heat transport equation to study the 
effects of radial heat flux on the dynamics of collapsing fluid sphere. From 
this coupled equation f p7|) . we flnd that the Gauss-Bonnet term decreases 
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(increases ) the value of a for which Fg^av has larger (smaller) values which 
increases the collapse of the system. Also, the equivalence principle is valid 
in /(G) gravity. We have explored that when a exceeds the critical value, the 
gravitational force plays the role of antigravity and reversal of collapse would 
occur in the system. A model with bouncing behavior has been presented 
numerically by Herrera et al.^^-* We would like to mention here that the 
homogeneity in energy density and conformal flatness for any metric are 
necessary and sufficient conditions for each other, when j[G) is constant 
and system under consideration is composed of isotropic perfect fluid. This 
analysis can be extended for generalized f{G) model^^-* like /(G) = C\G + 
6*2^4^, Gi, G2 and a are constants. 



Appendix 



The interior metric ( !T2|) has the following non-zero components of the Rie- 
mann tensor 

i?oioi = XX" + - ^x'y' + ^xr , 

, , Y XY fxVx'Y' 1, 



i?i2i2 = {Yr 



R 



2323 




-^0303 = sin 6R0202, -R0313 = sin 9Ro2i2, -R1313 = sin 6Ri2i2- (43) 



This implies that the Riemann tensor has five independent components. The 
Kretchman scalar becomes 

^ ^ ^ ((XF)4^^°^°^''^ ^ (XrF)4^^°^°^''^ ~ (Xr3r2)2^^02i2)^ 



(^2^)4(^1212)'+ ^^^^g ^.^2^- 



+ 7^7^(^1212)' + 7^73^— T7(^2323)M- (44) 
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The Riemann tensor components in terms of the Einstein tensor and mass 
function can be written as 

-Roioi = {XYY ( Coo — 7r77?Cii + , ■.r, G22 



2X2 2^2^- (rr)2 (Yr) 



(1 in 

-K0212 — — ^ — ^01, 



2X2 ^Yry 

^2323 = 2mrrsin2e. (45) 
Inserting these in the Kretchman scalar fj44l) . it follows that 
48m2 16m /Goo Gn G22 



n 



{Yrf (Yrf Y"^ {rY) 



Gm \ , n ( f Goo\ , f Gii\ \ . ^ f G22 



3 + ^] +4 



XY J IV^V J I \{rY) 

_ GoqGii . I Gqo ^ Gil \ Gil / 
X2F2 Vx^'Y^ylp^- ^ ^ 
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